FORC Analysis of homogeneous nucleation in the two-dimensional kinetic Ising model 
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The first-order reversal curve (FORC) method is applied to the two-dimensional kinetic Ising 
model. For the system size and magnetic field chosen, the system reverses by the homogeneous 
nucleation and growth of many droplets. This makes the dynamics of reversal nearly deterministic, in 
contrast to the strongly disordered systems previously studied by the FORC method. Consequently, 
the FORC diagrams appear different from those obtained in previous studies. The Kolmogorov- 
Johnson-Mehl-Avrami (KJMA) theory of phase transformation by nucleation and growth is applied 
to the system. Reasonable agreement with the Monte Carlo simulations is found, and the FORC 
method suggests how the KJMA theory could be extended. 

PACS numbers: 75.60.-d, 77.80.Dj, 64.60.Qb, 05.50.+q 



I. INTRODUCTION 

The study of hysteresis loops can produce significant 
dynamical information, but it reduces the often complex 
dynamics of magnetization reversal to only a few quan- 
tities, usually the hysteresis loop area and the coercive 
field. The first-order reversal curve (FORC) method was 
recently developed^ to extract more information from ex- 
periments on magnetic systems. It has produced inter- 
esting results, mostly in systems with strong disorder in 
the physics of magnetismSj&l and in the geosciencesi^ 



II. 



FORC ANALYSIS AND MONTE CARLO 
SIMULATION 



The FORC technique involves saturating the magne- 
tization M in a positive field Ho, decreasing the field to 
a series of progressively more negative return fields H r , 
and then in each case increasing the field back to Ho 
while recording the magnetization m = M/M s , where 
M s is the saturation magnetization. This results in a set 
of curves m(H r , H), where H represents the field as it is 
increased from H r back to Ho, as shown in Fig.^a). 

We use FORC analysis to better understand the pro- 
cess of hysteresis in the two-dimensional Ising model, 
which is defined by the Hamiltonian 



H = -J^2SiS s -nH(t)^2Si 

(i,3) ' 



(1) 



where the lattice spins Si — ±1, the exchange constant 
J is in units of energy, fi is the magnetic moment per 



site, and H(t) is the applied time-dependent field. The 
simulation is run on an L x L square lattice with peri- 
odic boundary conditions. The first sum is taken over 
all nearest-neighbor pairs of spins, and the second sum 
over all individual spins. Hereafter we use units such 
that J = /i = fcs = l. Our simulations are run at a 
temperature of T = 0.8T C , and we use the Glauber dy- 
namic with random, single-site updates for the Monte 
Carlo (MC) simulation^ The field decreases linearly at 
the sweep rate — Q from Ho to H r , and then increases 
linearly at the sweep rate +Q back to Ho- 

The free-energy barrier for the metastable state in the 
Ising system leads cither to a co-existence (CE), single- 
droplet (SD), multi-droplet (MD), or strong-field (SF) 
regime for magnetization reversal. 8 Occurring at inter- 
mediate values of system size, field and temperature, the 
MD regime is arguably the most accessible to experi- 
ments. The SF regime requires fields that may be too 
large to generate experimentally, while the SD regime re- 
quires small systems or very low temperatures, and thus 
very long metastable lifetimes. To ensure that the rever- 
sal is in the MD regime, the saturation field is chosen as 
Ho = 0.545, and the system size as L = 512. We ob- 
served directly that at L = 512 more than 100 droplets 
were involved in the reversal, and our FORC results did 
not change appreciably for larger L or larger Hq . We cal- 
culated the characteristic reversal time in a constant field 
H = —Ho, defined as the time to decrease from m = 1.0 
to m — -0.8, to be roughly 100 MCSS (MC steps per 
spin). We therefore chose fi = 2.18 x 10" 3 MCSS" 1 , 
so that the time for the field to decrease from H = 
to H = —Hq was long enough to produce full reversal. 
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In Fig. ^a) we present the resulting family of FORCs, 
where each FORC is an average over 20 MC realizations. 

It is useful to generate from this family of FORCs the 
FORC distribution, 

ld 2 m(H r ,H) 

^ H ^ = -2^wMr- (2) 

The motivation for this definition is that when the 
FORCs are taken from a Preisach system, the FORC dis- 
tribution reproduces the Preisach distribution^ In order 
to calculate the FORC distribution from our MC data, 
as in Ref. |jj we fit the family of curves m(H r , H) around 
each point by a second-order polynomial g(H ri H) = 
aH r H + bH 2 + cH 2 + dH r + eH + /, using a linear least- 
squares method. The FORC distribution is most clearly 
displayed in a contour plot of its level curves, called a 
'FORC diagram,' which we show in Fig. |2ta). 

There are several interesting features of the FORC re- 
sults. First, a sizable part of the FORC diagram is neg- 
ative. Since the FORC distribution measures how the 
slope dm/dH varies from one FORC to the next, which 
can be quite different for different dynamics, this is a rea- 
sonable result. It is only in the case of the Preisach model 
that the FORC distribution must by definition be posi- 
tive. Second, the positions of the minima of the FORCs, 
naively expected to be at H = 0, gradually shift from 
a small negative field to a small positive field. We will 
provide insight into these features in section ITTTl 

III. APPLICATION OF THE KJMA THEORY 

We briefly sketch the idea of the Kolmogorov-Johnson- 
Mehl-Avrami (KJMA) theory. See Refs.||and|! for more 
complete descriptions. In the present case of MD rever- 
sal, the system parameters are related as a <C r c <C ro <C 
L, where these parameters represent the lattice spacing, 
the critical droplet radius, the average distance before 
droplets meet during growth, and the system size, re- 
spectively. The critical droplet radius is the radius be- 
low which droplets of overturned spins tend strongly to 
shrink, and above which they tend strongly to grow. 

For the MD regime, magnetization reversal can be ap- 
proximated by the KJMA theory of droplet nucleation 
and growth, in which critical droplets are assumed to nu- 
cleate and grow with field- and temperature-dependent 
rates. First, the extended volume <E>(i) of the expanding 
droplets, uncorrected for droplet overlap, is calculated. 
Then the expression <f>(t) = exp[— $(£)] is applied to find 
the volume (j)(t) of overturned spins, corrected for droplet 
overlap. (See Ref. and references cited therein.) 

To apply the KJMA theory, we included the linearly 
varying field in the KJMA expressions, found straight- 
forward equations for the extended volume in the region 
of negative fields, and then numerically integrated these 
equations. As in Sides et al.j*» we approximated the be- 
havior of the system in positive fields as a combination of 
a shrinking negative metastable region, whose volume is 



estimated as the time reversal of the calculated growth in 
negative fields, and nucleation of a positive stable phase 
within this region, estimated from the calculated nucle- 
ation in negative fields. Because the calculated growth 
in negative fields extends only to H — H r , this approxi- 
mation can continue only up to H — \H r \. However, the 
system at H = \H r \, where positive droplets are grow- 
ing, can be approximated by a similar (but oppositely 
magnetized) state of the system, obtained from the full 
hysteresis loop. We found the point in the calculation of 
the full hysteresis loop where the magnetization is equal 
in magnitude to, but opposite in sign from, the magne- 
tization at H = \H r \. Beginning with the results of the 
KJMA integration at this point, we continued the inte- 
gration from H — \H r \ to H = Hq. 

In comparing the calculated results to the MC simula- 
tion data, we note first that the approximation of treating 
the final portion of each FORC loop using results from 
the KJMA simulation in negative fields works well, de- 
viating only in small kinks (almost invisible on the scale 
used here) which develop in the first FORCs in Fig. QJb) . 
This suggests that the behavior of the remaining pos- 
itive metastable phase in negative field values is simi- 
lar to the behavior of the remaining negative metastable 
phase in the last stage of the FORCs. The first FORCs 
are least similar, probably because they consist of sev- 
eral well-defined shrinking droplets, rather than an irreg- 
ularly spaced and shaped volume. It is also interesting 
that the magnetization decreases from saturation much 
sooner in the MC family of FORCs, due to the presence 
of subcritical droplets which are neglected in the KJMA 
treatment. This causes the values of the FORC distribu- 
tion to be larger for the KJMA case, since the slopes of 
the FORCs must be higher. 

We produced animations of the MC simulation in order 
to better understand the dynamics. In Fig. da), the 
FORCs beginning in the range —0.2 > H r > —0.25 have 
their minima at H < because droplets that are just 
above the critical radius become subcritical and shrink 
as the field increases back toward H = 0. The FORCs 
beginning in the range —0.25 > H r > —0.32 have their 
minima at values H > for a similar reason. These 
minima occur at m < 0, where the magnetization can 
be viewed as consisting of positive droplets in a negative 
background. In the subsequent small positive fields, the 
growth of the larger positive droplets is accompanied by 
the shrinking of the smaller, subcritical positive droplets. 
The entire shift of the minima is reflected in the presence 
of values p < of the FORC distribution in the region 
(-0.30 < H r < -0.25, -0.05 < H < 0.05) in Fig. 

In areas of the FORCs in which a full reversal occurred, 
there is a loss of memory in the MC simulation, resulting 
in a value p ps for these areas in Fig. |2Ia). We also 
studied the families of FORCs and FORC diagrams for 
the same system with ti = 2.18 x 10~ 4 MCSS" 1 (not 
shown) . We found the form of the MC and KJMA results 
were quite similar to the Q = 2.18 x 10~ 3 MCSS" 1 case, 
except that the coercive field was about half as large. 
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FIG. 1: Family of FORCs for L = 512 kinetic Ising model at temperature T — 0.8T C in a linearly varying field with sweep rate 
ft = 2.18 x 10~ 3 MCSS -1 . (a) Data from MC simulation, (b) Calculation from application of KJMA theory. 
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FIG. 2: FORC diagrams for the system described in Fig.0 The magnitudes in the contour plot range from p = —400 (black) 
to p = (large gray area) to p — +400 (white), (a) Calculated from MC simulation data, (b) Calculated with KJMA theory. 



IV. CONCLUSION 

In conclusion, we have applied FORC analysis to the 
two-dimensional kinetic Ising model, finding reasonable 
agreement between the predictions of KJMA theory and 
the results of MC simulation. The FORC analysis offers 
insight into the discrepancies between the KJMA and MC 
results, which we have identified as due to the effects of 
subcritical droplets and loss of memory. 

In future work, we hope to extend the KJMA theory 
to include these effects and to make predictions with the 
KJMA theory for longer time scales, where MC simula- 



tion is not practical. Although the simulation results for 
ft = 2.18 x KT 3 MCSS" 1 and ft = 2.18 x KT 4 MCSS" 1 
were very similar, the long time scales of physical systems 
may produce different and physically relevant insights 
into the dynamics of the homogeneous reversal process. 
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